Free and bound spin-polarized fermions in the fields of Aharonov— Bohm kind 



V.R. Khalilov, I.V. Mamsurov and Lee Ki Eun 

Faculty of Physics, Moscow State University, 119991, Moscow, Russia 

The scattering of electrons by an Aharonov-Bohm field is considered from the viewpoint of 
quantum-mechanical problem of constructing a self-adjoint Hamiltonian for the Pauli equation. 
The correct domain for the self-adjoint Hamiltonian, which takes into account explicitly the elec- 
tron spin is found. A one-parameter self-adjoint extension of the Hamiltonian for spin-polarized 
electrons in the Aharonov-Bohm field is selected. The correct domain of the self-adjoint Hamil- 
tonian can contain regular and singular (at the point r = 0) square-integrable functions on the 
half-line with measure rdr. We argue that the physical reason of the existence of singular functions 
is the additional attractive potential, which appear due to the interaction between the spin magnetic 
moment of fermion and Aharonov-Bohm magnetic field. The scattering amplitude and cross section 
are obtained for spin-polarized electrons scattered by the Aharonov-Bohm field. It is shown that 
in some range of the extension parameter there appears a bound state. Since the Hamiltonian of 
the nonrelativistic Dirac-Pauli equation for a massive neutral fermion with the anomalous magnetic 
moment (AMM) in the electric field of a linear charge aligned perpendicularly to the fermion motion 
has the form of the Hamiltonian for the Pauli equation in the Aharonov-Bohm flux tube, we also 
calculate the scattering amplitude and cross section for the neutral fermion. 
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I. INTRODUCTION 

The quantum Aharonov-Bohm (AB) effect, predicted in is an important physical phenomenon. 
Due to the cylindrical symmetry of Aharonov-Bohm field configuration, the relevant quantum mechanical 
system is invariant along the symmetry (z) axis and essentially two-dimensional in the xy plane. So, the 
models applied to describe AB effect can usually be reduced to the (2+l)-dimensional ones. Such models 
are studied by means of the Dirac equation in 2+1 dimensions and results obtained are applied to other 
problems. Solutions to the two-component Dirac equation in the AB potential were first obtained and 
discussed by Alford and Wilczek in Ref. 2] in a study of the interaction of cosmic strings with matter 
and the production of particle-antiparticle pairs by the vector potential of a moving cosmic string. The 
effect of vacuum polarization in the field of infinitesimally thin solenoid is recently investigated in [3j and 
a wonderful phenomenon is revealed: the induced current is finite periodical function of the magnetic 
flux. 

The scattering of spinless charged particle by an infinitely long thin solenoid was examined by Tyutin in 
Q from the viewpoint of quantum-mechanical problem of constructing of self-adjoint Hamiltonians. The 
Hamiltonian for the Aharonov-Bohm problem is essentially singular and hence it cannot be immediately 
defined in the domain [0, oo) for any differentiable and enough rapidly decreasing (at r — > oo) functions 
in the Hilbert space of square-integrable functions on the half-line with measure rdr [J. Usually, the 
Hamiltonians are symmetric operators in its natural domain. The problem of constructing a self-adjoint 
Hamiltonian is to find all self-adjoint extensions of given symmetric operator and then to select correct 
self-adjoint extension. The correctness of the known Aharonov-Bohm solutions for scattering problem 
was analyzed for spinless particles in Q , in which a self-adjoint extension of the Hamiltonian is selected 
by physical condition - "the principle of minimal singularity" : the Hilbert-space functions for which the 
Hamiltonian is defined must not be singular. 

A one-parameter self-adjoint extension of the Dirac Hamiltonian in 2+1 dimensions in the pure AB field 
was constructed by means of acceptable boundary conditions in [5|-|8|. For some extension parameters 
a domain of the self-adjoint Hamiltonian was found to exist in which the Hamiltonian is self-adjoint 
and commutes with the helicity operator One- and two-parameter families of self-adjoint Dirac 
Hamiltonians in the superposition of the Aharonov-Bohm solenoid field and a collinear uniform magnetic 
field in the respective 2+1 and 3+1 are constructed in (Tol. Hlj. 

Nevertheless, there remains a question about the correctness of solutions for the scattering of spin 
particles in the AB flux tube. Indeed, for this problem we need to use the Hamiltonian for the Pauli 
equation, which contains one more extremely singular (spin) term taking account of the interaction 
between the spin magnetic moment of electron and magnetic field. 

An interesting and important corollary to the Aharonov-Bohm geometric phase is a phase acquired 
by the wave function of a neutral massive fermion with the magnetic moment when it propagates in an 
electric field of a uniformly charged long conducting thread aligned perpendicularly to fermion motion 
(such a field will be called the AC configuration). The fermion transport is influenced by the phase 
acquired by the fermion wave function and the resultin g ph ase difference leads to a spin-field dependent 
effects in scattering (the Aharonov-Casher (AC) effect [lja]). The Hamiltonian of nonrelativistic Dirac- 
Pauli equation for a massive neutral fermion with the anomalous magnetic moment in the electric field 
of a linear charge aligned perpendicularly to the plane of fermion motion has the same form with the 
Hamiltonian for the Pauli equation in the pure AB potential. The nonrelativistic quantum motion of an 
uncharged massive fermion with an anomalous moment interacting with the electric and magnetic fields 
produced by linear electric and magnetic sources in a conical spacetime was studied in [13j. 

Thus, the self-adjoint extensions for the nonrelativistic Hamiltonians with the particle spin (this term 
leads to additional singular 6(r) potential in the Hamiltonian) have still not been constructed. 

This paper is organized as follows. In Section II we find a one-parameter self-adjoint extension of the 
Hamiltonian for spin-polarized electrons in the Aharonov-Bohm field. In Section III we find the scattering 
amplitude and cross section as well as the wave function of bound state for spin-polarized electrons in 
the Aharonov-Bohm flux tube of zero radius. In Section IV we discuss the scattering problem for spin- 
polarized neutral massive fermion with AMM in the field of uniformly charged conducting thread of zero 
radius. In Section V we briefly discuss physical results. 

We shall adopt the units where c = K = l. 
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II. SELF-ADJOINT EXTENSIONS FOR THE PAULI HAMILTONIAN IN AN AB 

POTENTIAL 



The Pauli equation for an electron of mass m and charge e < in an AB potential specified in cylindrical 
coordinates as 



A° = 0, A r = 0, A v = -, A z =0, B=^-, (1) 

r Z7T 



r = \J x 2 + y 2 , tp — arctan(y/cc) 

is 

i^(t,r,z)=HV(t,r,z), r = (x,y), (2) 
where ^(t, r, z) is a spinor and the Hamiltonian H is 



Here 03 is the Pauli spin matrix 

0-3 



1 

-f 



The last term in describes the interaction of the electron spin with a magnetic field. 

Equation (fTJ) is the potential of an infinitely thin solenoid with a finite magnetic flux $ in the z direction 
in the range r > 0. The magnetic field H is restricted to a flux tube of zero radius 

H = (0, 0, H) = V x A = Bw5(r). (4) 

Therefore, we need to construct a self-adjoint Hamiltonian of Eq. in cylindrical coordinates for the 
electron motion in the xy plane. It is seen that the inclusion of particle spin in the AB problem leads to 
the additional 6(r) potential in the Hamiltonian [14| . 

We seek solutions to the Pauli equation for the electron motion in the xy plane in the form 

oo 

*(t, r, <p) = exp(-iEt)f(r, <p)ij> = exp(-iEt) V" Fi(r) e-x^(ilip)ip s , (5) 

I — — OO 

where E is the electron energy, I is an integer, and ip s is a constant two-spinor 



^~ J 2 \ 1 - a 

where s = ±1 characterizes the electron spin projection on the z axis. 
The radial part Fi (r) of the wave function satisfies 

d 2 Id I 2 2\e\Bl e 2 B 2 nri .„ S(r)\ , . 

-Ti+—, 2-— 2 j — 1~ 2Em - efls— ) Fi(r) = (7) 

and it depends explicitly on the number s, which selects a particular value of the spin projection on the 
z axis. 

At first, it is rewarding to remind the problem of constructing solutions for the linear differential Bessel 
expression for functions Fi(r) or f(r) = F[(r)/y/r of the forms 

\ dr z r ar r z J 

or 

where /i = |e|B. In the form ([9|) the quasiderivative 1(f) is self-adjoint and can be defined in the 

range — oo < x < oo, x = r. 



4 



In the Hilbert space £2(0, b), i.e. the space square integrable functions in the range (a, b) (the interval 
(a, b) may be infinite) the linear everywhere dense set D' is considered. The set D' consists of the infinitely 
differential functions that vanish outside some cut [c, d] (the cut may be various for various functions) 
contained completely in the range (a, b). On D' Q the operator L is defined by equality L f = 1(f) (l5lll6j. 
In each cut [c, d] the Lagrange identity 



d d , 

J Kf)9dr - J fl(g)dr = (/f - |«) = [/,*]< 



(10) 



is valid for any function g being complex conjugate of Q-So Lq is a symmetric operator. The closure Aq 
of operator Lq will be also a symmetric operator [ID, [l6j]. We need to construct self-adjoint extensions 
A of operator Aq. The adjoint operator Aq is given by equation 

A* f = l(f). (11) 

is defined on the set D(Aq) of all functions that are absolutely continuous with their first derivative 
in the range (a, b) and the quasiderivative 1(f) belonging to L 2 (a, b). 

If Eqs. (|8]), (O are considered in the range (0, 00) then they will be singular on both ends. In this case 
the domain D(A ) of operator A contains all functions / from D(Aq) for which 

[f,9t = (12) 

for all g from D(Aq). In the most important case fj, is nonintegral and the corresponding operator A 
will have the deficiency indices (2,2) at < \l + fx\ < 1. The eigenfunctions G (r) = g^y/r of the adjoint 
operator Aq with the eigenvalues ±i satisfy equations 

-^ + ^-^T^Wm = 0, (13) 
dr z r ar r- 2 / 

where p 2 7^ is inserted for dimensional reasons. Square integrable normalized solutions of (|13[) are 
expressed by the Bessel functions of third kind (see, for instance (l7l fl8j) 

G+(r) = A+H^^^pr), G7(r) = A; (e^^pr), (14) 



where Af 1 are the normalization factors and 



u = \l + fj,\. (15) 



If fi is nonintegral then the operator Aq has two eigenfunctions belonging to £2(0, 00) for each eigenvalue 
±i. Indeed, writing fi as /x = \p] + 7, where [/x] is the largest integer < /x, it is seen that these functions 
belong to the critical subspace < 7 < 1 (with I — — [[i] — 1). The self-adjoint extensions of Aq can be 
parameterized by the parameter 8. The correct domain D(A ) for the self-adjoint extension A 8 of Ao is 
then given by 

D(A ) = D(Aq) + dif^e^V) + C 2 H^(e^' 4 pr) + Bi^^e^V) + B 2 H^ (e^^pr), (16) 

where Bk = Cfe exp(— iQ), Ck are the arbitrary complex numbers and 6 is an arbitrary (but fixed for a 
given extension) parameter. The domain D(Aq) contains absolutely continuous, square integrable on the 
half-line with measure rdr and regular (at the point r = 0) functions. If /i is a non-integer then the 
operator Aq has two eigenfunctions belonging to £2(0, 00) for each eigenvalue ±i. Let us write (see, fl8j ] ) 

mx) = j[e-^Mz)-j.^)] t mz)ss _i[*"Mz)-J-A*)] t (17) 

sin vtt sin vtt 

where J v (z) and J_ v (z) denote the regular and irregular Bessel functions. The asymptotic behavior of 
the Bessel function at r — > is 

JJx) -, (18) 

where T(z) is the gamma function of argument z. 
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Applying (fTS)) we derive 



CHjPieWpr) + Ce- l6 Hi 2 \e^/ 4 pr) = ^ x 

sm -KV 

in which 1/ = 7, 1 — 7. 

For any function f(r) from the defect subspaces the Lagrange identity (|10[) has to be satisfied 



(19) 



(20) 



Since the functions from defect subspaces must be square integrable on the half- line [0, 00) with measure 
rdr (i.e. they, in fact, obey the Lagrange identity at r — > 00) it is enough to put a boundary condition at 
the origin 



Hi,, I f^l _ tf_g 

r->o \ dr dr 



0. 



(21) 



A symmetric operator is self-adjoint if its domain coincides with that of its adjoint (see, for instance 
Q). So, one has to posit the same boundary condition at the origin for functions from the spaces G ± (r) 
defined on the half-line. 

The correct domain for the self-adjoint extension of the Hamiltonian of spinless particles was selected in 
[H by the physical condition of "minimal singularity" - the domain of the Hamiltonian should be composed 
by entirely regular functions. It follows from (fT9|) that this requirement fixes the parameter 9 = irj/2 
for v = 7 and 9 = ir(l — 7)/2 for v = 1 — 7. Since < 7 < 1 the correct domain of the Hamiltonian is 
composed by entirely regular functions for < 9 < n/2. Insisting on regularity of all functions at the 
origin forces one to reject irregular (but square integrable on the measure rdr) solutions in the defect 
subspaces, which can entail a loss of completeness in the angular basis. Indeed, the inclusion of particle 
spin in the AB problem leads to the presence of spin-dependent potential term (~ (\e\sB/r)5(r)) in the 
Hamiltonian, which can be taken into account also by means of the boundary condition at the origin. Let 
B be chosen, for definitencss, positive. Then this potential is repulsive for s = 1, and attractive for s = — 1 
what means that if s — 1 then functions from the defect subspaces defined on the real half-line should 
be regular at the point r = but if s = — 1 then singular (concentrated at the origin) square integrable 
functions have to appear. It should be emphasized that the appearance of concentrated solution at the 
point r = is due to the physical reason, namely, by the presence of attractive potential (|e| sB/r)S(r). 
We see from (fl~9|) that the defect subspaces does not contain the regular functions with indices 7 and 
1 — 7, respectively, at 9 = 3wy/2 for v = 7 and 9 = 37r(l — r y)/2 for v = 1 — 7. The range of the extension 
parameter in this case is tt/2 < 9 < 3tt/2. 

It should be noted that the linearly- independent solutions of equation 



dr 2 



fir) = \f 



with the Bessel operator (j9]) are expressed as 



f(r) = aiy/rJ u (V\r) + biy/rJ- v {V\r), 



(22) 



(23) 



where a/ and bi are the normalization factors and A a is a real positive number. The continuous part of 
spectrum of any self-adjoint extension of operator Aq (|22p considered in [0,oo) coincides with the positive 
semiaxis A > 0. 

The self-adjoint extensions of the operator Aq are parameterized by (fl"6)) . For s — 1 the domain D(Aq) 
is composed of regular (for example, the regular Bessel J^vAr)) functions for the range of the extension 
parameter < 9 < n/2. Because of the invariance of Eq. ([7} under simultaneous transformations 
(i — > — fi, I — > — I and s — > — s we shall restrict ourselves to study of the case /1 > 0. 

Therefore, for s — 1 the complete set of functions of the self-adjoint Hamiltonian is given by the regular 
functions defined on the real half-line 



f(r,<P) = J\i +ft \(V\r)exp(il<p), I = 0,±1,±2,.. 



(24) 



For s — — 1 for the range of the extension parameter tt/2 < 9 < 3ir/2we derive that the defect subspaces 
contain two singular Bessel functions J_ 7 (\/Ar) and J 7 _i(v / Ar) defined on the real half-line but bearing 



G 



in mind that < 7 < 1 we have to leave only one of them. We shall leave J_ 7 (vAr). In order to the set 
of functions would be complete in I, we need to remove from the domain D{Aq) the function J 7 (vAr) 
with I = — \p]. Hence, the complete set of functions of the self-adjoint Hamiltonian is given by 

f(r, tp) = J, ;+Al |(VAr) exp(il<p), I = 0, ±1, ±2, 1 ^ -[/x], . . . ; J^(V\r) exp(-i\ji](p). (25) 

Therefore, the complete set of functions of a self-adjoint Hamiltonian for the Pauli equation in the AB 
potential is a set (j24|) for s — 1 and a set (|25|) for s = — 1, in which wc must replace in the argument of 
Bessel functions V\ by k = \j2mE. 

In addition, if the self-adjoint extension of Aq is selected by the range of the extension parameter 
7r/2 < 9 < 3tt/2 there appears a bound state, which is expressed through the MacDonald function 
if(i_ 7 )(z) (or Kj(z)). The appearance of bound state "suggests that this range of parameters in the 
effective Hamiltonian parameterizes nontrivial physics in the core" @. One sees from our model the 
appearance of bound state is possible only if the additional spin term in Eq. ([7]) is attractive, i.e. s = — 1 
(see, also [3 (3). 



III. SCATTERING OF ELECTRONS BY AN AB FIELD 



We shall consider in what follows the case 

/i=[ M ]+7 = n + 7) (26) 

where n is an integer and 

1 > 7 > 0. (27) 

The expansion of the spatial electron wave function is given by 

00 

il>{r,(p)= J2 NiJ v {kry 1 ^, s = l, (28) 

/— — OO 

where Ni are constants and 

00 

i/)(r,<p)= J2 NiMkry 1 * + A=-„J- 7 (fcr)e-»^, s = -1. (29) 
/— — 00 

Here the summation is carried out with the omission of the I = —n term. 

If we assume that the electron wave propagates from the left along the x axis so as to incident wave 
is tp = e lkx , then ip is the scattering angle measured from the positive x axis. The asymptotic form 
of electron wave function at r — > 00 is a superposition of ingoing plane wave and scattered outgoing 
cylindrical wave 

i> k (r, if) = e lkx + lM e i(kr-n/4) (30) 
At r — > 00 the plane wave is given by 

gikreosip _^ \^ e il<P |gi(fer-7r/4) _|_ (_]\l e -i{kr— 7r/4) \ 

\j2ixkr , \ < 

I—— 00 

and Eq. ([28)) is given with using the expansion of Bessel functions as follows: 
1 00 

TP( r>(p ) y Nie U( V +*/2) ( e i{kr-*/A) e -m\l+»\/2 + e -i(kr-ir /4) e ^\l+^\/2\ (33) 

y/2irkr - ' 



/— — 00 



Choosing the coefficients Ni so as to the ingoing waves would be canceled in Eq. ([30]). one easily obtains 

Nl = e -i(*/2)|J+/i|. (33) 

The scattering amplitude for solution ([28]) is then found to be given by the AB formula 

e -i(n+l/2) V 

h{ip) = sin ^7 : — — . (34) 

y/znk sm((p/2) 
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Inserting in (j2"9")l the (±) exp(— itt(^ — n)/2)J 1 {kr) terms we obtain the scattering amplitude in the AB 
form again 

e -i(n-l/2)(p 

f-i(<p) = sin^ 7 r— . — — ■ (35) 
v2nk sm(p/2) 

Therefore, if electrons are polarized so as to their spins are oriented along the z axis the scattering 
cross section in an AB potential is given by the AB formula 

da_ = sin 2 tt7 
dp 2vrfcsin 2 (( / 9/2)' 

The cross section for unpolarized electrons is obviously described Eq. (|36[) too. 

It is easily to see that when the electron spin is perpendicular to the z axis we must take the wave 
function in the form 

oo 

^(r, ip) = NiJ v (kr)e ilv + e' mv (J 7 {kr)e- l7r ^- n ^ 2 + J_ 7 (fcr)e"^-™ )/2 ) /2, (37) 

/ — — oo 

where the summation is carried out with the omission of the I = —n term. 
Now the scattering amplitude has the form 

sin7P7 /e _4 (™ +1 / 2 ^ - \ 
Wy y^ffc V sin((^/2) y 1 ' 

and the cross section is 

da sin 2 7T7 



dip 2irk \sm 2 (<p/2) 



- 1 (39) 



in case the electron spin lies in the plane of scattering. In case | | < 1 this formula describes the scattering 
cross section of spin-polarized electrons by a long cylindrical magnetic flux tube of small radius 



A. On a bound electron state in the AB held 



It was shown in Section II that in the range of the extension parameter tt/2 < 9 < 3tt/2 there appears 
a bound state, which is expressed through the MacDonald function Kn_ 7 -j(z) (or Kj(z)). The radial 
part of wave function of bound state is solution of Eq. with s = — 1 and E < 0. It is Kn_^\(z) (or 

K~ ( (z)) of argument z = r^j2m\E\ dchncd on the real half-line where m is the electron mass and E < is 
the energy of bound state. One sees that a bound electron state may occur in the quantum system under 
consideration if the interaction of the electron spin with a magnetic field is included and it is attractive. 



IV. SCATTERING OF A MASSIVE NEUTRAL FERMION IN THE AC FIELD 

CONFIGURATION 

If we put |e| = Ms in equation J7| it will describe the motion of a massive neutral fermion with the 
anomalous magnetic moment M in the electric field of a linear charge aligned perpendicularly to the 
fermion motion. Then, v is 



v=\l + s~{\^Q, 7 = MB, k = y/2mE, (40) 

where B/2 is the linear charge density, s — ±1 is the fermion spin projection on the z axis, E is the 
fermion energy and the results obtained above are valid for a neutral fermion with AMM in the AC 
configuration. 

The scattering amplitudes for/z>0, /i<0, s = ±1 are respectively [2(| 

^ = 72^ ' (41) 



S e i»(\n\-l/2)<p+i\n\n gW^UI) 

/sM = "71^ ^72) ' (42) 

Thus, if the spin of neutral fermion is oriented along the z axis the scattering cross section in an the AC 
configuration is given by 

da = sin 2 (7T7) ^ 43 ^ 
ac 27rfcsin 2 i^/2 

The cross section for unpolarized neutral fermions with AMM in the AC configuration is also described 
by Eq. (@3j) Eqs. gl]), ((42j) and (j43]) coincide respectively with Eqs. d34j), §B and ([36) . 

If the spin of neutral fermion in the initial state lies in the scattering plane then the scattering amplitude 
and cross section are respectively 

M<P) = (~l) n+1 /o sin(n+l/2)y, (44) 

yzirk sm(ip/2) 

and 

^ = |/oM| 2 = -f^^sin> + l/2V. (45) 
dip 2nksm (ip/2) 

One can show that these formulas are described the case of negative fi too, if n replaced by \n\ — 1 + 0(fi), 
where 6(x) = [1 + sign(a;)]/2. It is seen if the particle spin in the initial state lies in the scattering plane 
the scattering amplitudes and cross sections for AB and AC effects are different. 



V. DISCUSSION 



If the dependence of cross section is to be studied upon the spin polarization of particles in the initial 
state, the detector should detect only those scattered particles for which the angles between the spin and 
momentum vectors are equal in initial and final states. In the AB configuration the magnetic field H is 
localized inside infinitely thin solenoid, so the Hamiltonian ([3]) does not contain the vector cr matrix in 
the range r > and the electron spin s oriented along a unit vector n is conserved. Since, the scattering 
occurs in the xy plane the spin vector of electron in initial state can effect the scattering amplitude and 
cross section only when it has a nonzero projection on the plane of scattering. So, in fact, the dependence 
on the (initial) spin of electron in the cross section determines and arises because the direction of motion 
of electron is changed as a result of scattering (see [l9[). 

In case a fermion with AMM moves in the xy plane, the Hamiltonian contains only the 03 matrix. 
Therefore, only the fermion spin projection on the z axis s is conserved. Like the AB case the spin vector 
of fermion in initial state can effect the scattering amplitude when it has a nonzero projection on the 
scattering plane but in the AC case the dependence on the spin vector in the initial state in the cross 
section determines not only by the change in momentum of fermion but also by the change in the spin 
orientation of scattered fermion (the quasiclassical precession of fermion spin) . 
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